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Abstract 

The problem of the microscopic description of excited states of the even-even open- 
shell atomic nuclei is considered. A model is formulated which allows one to go beyond 
the quasiparticle random phase approximation. The physical content of the model is 
determined by the quasiparticle time blocking approximation (QTBA) which enables 
one to include contributions of the two-quasiparticle and the two-phonon configura- 
tions, while excluding (blocking) more complicated intermediate states. In addition, 
the QTBA ensures consistent treatment of ground state correlations in the Fermi sys- 
tems with pairing. The model is based on the generalized Green function formalism 
(GGFF) in which the normal and the anomalous Green functions are treated in a 
unified way in terms of the components of generalized Green functions in a doubled 
space. Modification of the GGFF is considered in the case when the many-body nuclear 
Hamiltonian contains two-, three-, and other many-particle effective forces. 



PACS numbers: 21.60.-n, 74.20.-z 



1 INTRODUCTION 



One of the most widely used approaches which are applied to the description of excitations 
of the even-even atomic nuclei is the random phase approximation (RPA, see, e.g., Ref. pQ). 
Within this approximation the nuclear excitations are treated as the one-phonon states 
which are superpositions of the one-particle-one-hole (lplh) configurations. However, this 
approach is applicable only in the case when pairing correlations are not essential, i.e., 
strictly speaking, only for magic nuclei. Generalization of the RPA taking into account 
pairing correlations explicitly is the quasiparticle RPA (QRPA), where the excited states 
(phonons) are expanded in the two-quasiparticle (2q) configurations. Thereby the QRPA 
extends the range of the RPA to the open-shell (non-magic) nuclei. Nevertheless, despite the 
significant success of both the RPA and the QRPA, there are several reasons for development 
of the models going beyond these approximations. 

First of all, description of the nuclear excitations in terms of the one-phonon wave 
functions is justified only for low- lying states. At higher excitation energies, fragmentation 
of the one-phonon states becomes important. This means that in addition to the lplh or 
2q configurations the more complex configurations should be incorporated (see Ref. [2]). 
The role of the effects related to the complex configurations is well manifested, for example, 
in the theory of giant multipole resonances (GMRs). It is well known (see, e.g., Ref. |3j) 
that RPA and QRPA enable one to describe the centroid energies and total strengths of the 
GMRs. However, both models fail to reproduce the total widths of the resonances and their 
fine structure. The reason is that these characteristics of the GMRs are significantly affected 
by the complex (mainly 2p2h or 4q) configurations which form the spreading width of the 
resonance. 

Another direction of developing a nuclear structure theory is associated with the 
models in which ground state correlations (GSC) beyond the RPA and QRPA GSC are taken 
into account (see Refs. [3J EJ EI). It has been shown that the GSC caused by complex 
configurations play an important role in the theoretical description of the experimental data. 
In what follows we will refer to this type of GSC as the GSC2 in order to distinguish them 
from the GSC1 included in the RPA and QRPA. In some cases, the GSC2 can strongly affect 
the transition strengths and can even lead to the appearance of new transitions which are 
absent in the calculations including complex configurations in the excited states only, i.e. in 
the calculations neglecting this type of GSC and using a restricted basis (see Ref. |B|). 
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A variety of models have been developed to study the effects of complex configura- 
tions on the structure of excited states of the even-even atomic nuclei (in addition to the 
aforementioned papers see also Refs. El EH UH EH El EH] and references therein). Never- 
theless, until recently (see Refs. [1211131), the quasiparticle-phonon model (QPM, Ref. 0) 
developed by Soloviev and co-workers was the only working approach which consistently 
treats the complex configurations and the pairing correlations on an equal footing. It is 
no surprising that comprehensive studies of the excitations of the open-shell nuclei taking 
into account complex (mainly two-phonon) configurations at the microscopic level have been 
carried out only within the QPM. In view of this the development of other approaches in 
this direction is particularly important. 

The principal goal of the present paper is to generalize the model of Ref. jTOj by 
including the pairing correlations. This model was developed to describe the excited states of 
the even-even doubly magic nuclei taking into account 2p2h (more precisely, lplhg>phonon) 
configurations. The model is based on the Green function (GF) formalism. The GSC2 were 
completely included within the model approximations. In the framework of this model the 
calculations of the GMRs in magic stable and unstable nuclei have been performed. Some 
of the results are presented in Refs. [31 EI] • I n these calculations reasonable agreement with 
the experimental data for the integral characteristics of GMRs, including the total resonance 
widths, has been obtained. Thus, one can expect that the extension of the model ^Oj to the 
open-shell nuclei will also give reasonable results. 

The second goal of the paper is to provide a modification of the GF formalism for the 
Fermi systems with pairing which is most suitable for solving the problem under considera- 
tion. To this aim, the generalized Green function formalism (GGFF) is presented, in which 
the normal and the anomalous GFs are treated in a unified way in terms of the components 
of generalized GFs in a doubled space. Modification of the GGFF is considered in the case 
when the many-body nuclear Hamiltonian contains two-, three-, and other many-particle 
effective forces. 

The paper is divided into two main parts. The first part (Sec. EJ) reviews the basic 
formulas and equations of the GGFF. The single-quasiparticle basis functions, which provide 
suitable representation of the model equations, are introduced. The second part (Sec. EJ) 
contains the formulation of the model in which pairing correlations, 2q, 2q(g>phonon, and 
two-phonon configurations are included. The model is analyzed within sum rule approach. 
The conclusions are given in the last section. 
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2 GENERALIZED GREEN FUNCTION FORMALISM 



2.1 Basic definitions 

Let a'(x) and a(x) be creation and annihilation operators of particles (free fermions) in 
the coordinate representation of the usual single-particle space. Here symbol x = {r, a, r} 
includes the spatial coordinate r, the spin a, and the isospin r variables. Considering the 
Fermi systems with pairing correlations it is convenient to pass from this single-particle space 
spanned by the coordinates x to the extended (doubled) space spanned by the coordinates 
y = {x,x}, where \ = ±1 is an additional index introduced for denoting the different 
components of the single-particle functions in the extended space (see Refs. PIE] for details). 
Let us define the operators b(y) = b(x, x) by the relations 

b(x,+) = a(x), b(x,-)=a\x). (2.1) 

From this it follows that b\y) = b(y) where y = {x, — x}- The Heisenberg representation of 
the 6-operators (in units where Planck's constant h = 1) reads: 

$( z ) = e lHt b(y) e~ lHt . (2.2) 

Here and in the following z = {t,y}, t is the time variable, H is a many-body Hamiltonian 
of an interacting fermion system. Obviously, these ^/-operators possess the property: 

&(z) = *{z), (2.3) 

where z = {t, y}. 

We will assume that the motion in the fermion system is determined by the nonrela- 
tivistic Hamiltonian H of the form 

H = H° + V, (2.4) 
where H° is a single-particle Hamiltonian including the external anomalous pair potentials: 
H° = / dx 1 dx' x \ h°(x 1 , x[) (x x ) a(x[) 



+ \ A°(x 1 ,x[)a t (x 1 )a f (x' 1 ) - \ A x[) a(xj a(x[) j , (2.5) 
V is an interaction including two-, three-, and other many-particle effective forces: 

K 

V = Y,V {k) , (2.6) 

k=2 
4 



— — J dx l ■ . . . ■ dx k dx 1 ■ . . . ■ dx k v *• \x 1 , . . . , x k ; x x , . . . , x k ) 

x a)(xi) ■ . . . ■ a\x k ) a(x' k ) ■ . . . ■ a(x[) . (2.7) 

Here and in the following J dx means the space integral over r and the sum over a and r 
indices. Analogously, in the following J dy will denote J dx and the sum over \i J dz will 
denote J dt dy. In case of the exact nuclear Hamiltonian we have: 

h°(x,x') = -(^i + fi T ^6(x,x'), A°(x,x') = 0, (2.8) 

where S(x,x') = 5{r — r') S a>a i S T y, \i T is the chemical potential for the nucleons with the 
isospin projection r which is introduced to simplify the following equations. Notice that the 
Hamiltonian H° can be formally rewritten in terms of the 6-operators as 

H° = \ J dydy"H°(y,y') b\y)b{y') + e , (2.9) 



where 

H°(x, + , x', +) = h°(x, x') , H°(x, + , x', -) = A°(x, x') , 

H°(x, - , x', +) = -A * (x, x') , H°(x, - , x', -) = -h°* (x, x') 



;2.io) 



e = - / dxh°(x,x) . (2.11) 
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Let | ) be the wave function of the ground state of the interacting fermion system. 
If the Hamiltonian H does not contain the external anomalous pair potentials A [e.g., if 
Eqs. f!2.8|) are fulfilled] , the number of particles is conserved, and | ) is an eigenfunction of 
the particle-number operator. However, in general case we shall not suppose that A = in 
Eq. ()2.5j) . i.e. we shall not suppose that the condition of the particle- number conservation 
is fulfilled for | ). Let us define the fc-particle generalized Green function (GF) in the time 
representation by the formula: 

G^(z l ,...,z k ;z[,...,z' k )=z^(0 | Tttfc) • . . . • *(z h ) ^(4) • . . . • ¥(z[) | ) , (2.12) 

where T is the time-ordering operator. In particular, for the single-particle GF we have: 

G{z, z') = G {1 \z ■ z') = -% ( |T*(z) ^\z') | ) . (2.13) 

The property 

G(z,z') = -G(z',z) (2.14) 

follows from Eqs. ()2.3|) and ()2.13|) . It can be seen from the definitions (|2.1j) and ()2.2|) 
that the normal and the anomalous GFs are the components of the generalized GFs G^ 
corresponding to the different values of the x-indices. 
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2.2 Equations of motion for the Green functions 



In case of the Fermi systems with pairing, the equations of motion for the many-particle GFs 
can be obtained with the help of the same technique based on the generating hmctionals 
depending on the auxiliary source fields which is frequently used for the Fermi systems 
without pairing correlations (see, e.g., Ref jH]). Let us define the generating functional W 
depending on the source field £ as 

W[£] =ln(0|TU|0) , (2.15) 

where 

U = «p( i /^ { (^»* W . M ). (2,6) 
It follows from Eq. (|2.3|) that one can consider the equality 

t(z,3>) = z) (2.17) 

to be fulfilled. Let us introduce the GFs with a source field £: 

r w (z , zM _ ( iTUggi) ■ • • • ■ *(* fc ) * f (4) • • • • • *M> I ) , 218 n 

lt^ \^z v . . . , z k , z v . . . , z k ) — i (0 1 TU 1 0) ' °> 

In particular, we have: 

/) , ; ,0 = -i (0|T ; V ;W )|0> = *) . (2.19) 

Obviously, coincides with G {k) defined by Eq. (jSH2D at £ = 0. 

It is easy to see that the GFs G^ can be obtained from the generating functional 
W[£] by a successive differentiation with respect to £. In particular, we obtain: 

8W 

G * M= w^y (2 - 20) 

L(2 z . z o _ jgj(£2 1 ji) 
where is the response function defined as 

L^(zi, z 2 ; z 3 , Zi) = Gf'(z 2 , z 3 ; zi, z A ) - G((z 2 , Zi) G^(z 3 , z 4 ) (2.22) 

(notice that this formula differs from the definition in Ref ^H] by the permutation of the 
arguments of L^). 



The equation of motion for the single-particle GF is obtained by the differentiation of 
G^(z\, z 2 ) with respect to time in analogy to the case of the Fermi systems without pairing 
correlations (see Ref EHI). It has the form 



(G°) ^i, z 2 ) + z 2 ) - £(z 2 , zi) = 1 (zi, z 2 ) + z 2 ) 



(2.23) 



where 



(G°) \z u z 2 )=[i 6(y 1} y 2 ) J- - H°(yi, y 2 ) ) Sfa - t 2 ) 



HlJy v') = $x,x' ( K X > x ')> ^£ * s ^ ne mass operator which is defined by the equations: 



dz"E^z,z")G^z",z') 



dz' 2 ■ . . . ■ dz' k dz'l 



W ^ (zi, . . . , z k ; z[, . . . , z' k ) 



A i l - k 
^ k\ (k- 1) 

k=2 V ' 

x G^ \zi , . . . , z k ; 2 , 2 2j • • • > ^/c) • 

s , b{t x - t; - xi ■ o) 

Al> Al 



X 



j=i 



M' X-l+l M> M + l 



X [$ x + i w ^( x l) ■ ■ ■ 1 X k ) X l) ■ ■ ■ 



5 X fc 



(2.24) 



(2.25) 



)]■ (2-26) 



In the Eq. (|2.26|) . is the antisymmetrized matrix element of the fc-particle interaction 
in the coordinate representation which is defined through the effective forces entering 
Eq. (12. 7|) and through the generalized antisymmetrized delta functions by the formulas: 



I; • • • j x k i -^1) • • • ) x fc 



X i 



k ! (k 



: ~ / c/x-i ■ . . . ■ dec i, dec-t 
- 2) ! J 



dx i. 



ry> ' rp' ' \ / ry>' ' ' iy' 



x \i ■ ■ ■ i x k 



X 



It ■ ■ ■ t X k 



u \x 1 , . . . , x k , x 1 , . . . , x k ) , 



(2.27) 



where 



£ | x ii ■ ■ ■ ? x k 
x ii ■ ■ ■ 5 X k 



( 8{x^ x' x ) . . . 5(x v x' k ) \ 



det 



V S(x k , x\) ... 8(x k , x' k ) J 
Notice that from Eqs. (l2~TUj) . (l2~T4l . (I2~m and (I2~231) it follows: 



(2.28) 



(2.29) 



The symmetry property of follows from its definition (|2.26j) : 

w 4, 4) = (-i)'w(*fe.., 4 ;%,■■■,**)■ (2.30) 

In order to obtain equations for the other (many-particle) GFs let us perform a 
change of the functional variable £ to the and consider a Legendre transformation of the 
functional W[£]: 

r[ G 5 ] = 2 iy [£] - j d Zl dz 2 [Z(z2, *i) - £(%, z 2 ) ] Gfa, z 2 ) . (2.31) 



Using Eqs. (gUj) and (j2~271|) we obtain 

5T 



^{z 1 ,z 2 )-^2,z 1 ). (2.32) 



The notation 5_ means that the variations of the GF conserve the property of antisymme- 
try (J2.14j) . This condition should be taken into account since variations of £, which generate 
the variations of G^, obviously do not lead to the violation of Eq. (j2.14J) . Conservation of the 
property (j2.14|) in the variational procedure can be automatically ensured if the following 
substitution is performed in a G^ -dependent functional: 

Gfc, z') = \[G^z, z') - Gs(z', z) ] . (2.33) 

In case of the vanishing source field the Eq. ()2.32|) leads to the stationarity condition: 

Using, further, Eq. ()2.23|) we obtain from Eq. ()2.32|) the following relation 

' — T = \ [ 1 ( z 2, Zi) G^ 1 (z 3 , z\) — G^ 1 (z~i, z 4 ) G^ 1 (z 3 , z 2 ) ] 

35 

5- G € ( zi , z 2 ) 5- G^ ( z 4 , z 3 ) ' 

Let us introduce an amplitude of the effective interaction 1^ which includes irreducible 
amplitudes both in the particle-hole (ph), and in the particle-particle (pp) channels: 

f(z 2 , Z\) 

H\Z\, Z2\ Z 3 , Zi) = i — — r. (2.36) 

0-Lr£{Zi, Z 3 ) 

From Eqs. (gHH) and (J2jgH) we obtain: 

Tt(zi, z 2 ; z 3 , z 4 ) = -2"e(z2, z\\ z 3 , z 4 ) = ~^(zi, z 2 \ z 4 , z 3 ) . (2.37) 



S 



In addition from Eq. (|2.35|) it follows that 

3-£\ z i, z l\ z 3, z i) — 2^(^4j z 3] z 2, z i) ■ (2.38) 
Notice that the response function defined by Eq. (j2.21|) satisfies the analogous equalities 

L&i, z 2 \ z 3 , z 4 ) = -L^(z 2 , Z\\ z 3 , z 4 ) = -L&i, z 2 \ z 4 , z 3 ) , (2.39) 

L^{z u z 2 ; z 3 , z 4 ) = L^(z 4 , z 3 ; z 2 , z x ) . (2.40) 

Differentiating Eq. ()2.23j) with respect to £ and then using Eqs. (j2.21|) and ()2.36j) . we 
obtain the Bethe-Salpeter equation (BSE) for the response function: 

L((zi, z 2 ; z 3 , z 4 ) = G^(z 4 , z x ) G^(z 2 , z 3 ) — G^(z 3 , z x ) G^(z 2 , z 4 ) 

- i J dz 5 dz 6 dz 7 dz$ G^(z 5 , z x ) G^(z 2 , z 6 ) 

X 2^5, z %\ z 7, Zb) ^O, Z ^ Z 3, z i) ■ ( 2 - 41 ) 

The equations for the many-particle GFs G^ with k > 2 are obtained by a differen- 
tiation of Eq. ()2.18|) with respect to £. Taking into account relation: 

— - = / dz 3 dz 4 L i {z x ,z% z 3 ,z 4 )—— r, (2.42) 

o£(zi,Z2) J 6Gi:{z 3l z 4 ) 

which follows from Eq. (|2.21|1 . we come to the recurrence formula 
Gc \z ± , . . . , z k ; z l7 . . . , z k ) = 



G^(z 1 ,z[) + J dz dz' L^(z[, z x ] z' , z) 

G$ 1 { z 2i ■ ■ ■ j z k i 4' • • • ) z 'k) ■ (2.43) 



5 

x 



Notice that in Eqs. ()2.42|) and (j2.43|) we have: 5_G^ = 5G^ owing to presence of the response 
function satisfying Eqs. (J2.39j) . 

The Eqs. (J222J), fl2H3), fl23U), and (f2~IHjl form the closed system of the 

functional differential equations of the GGFF. An important feature of these equations is 
that they do not change their form when the many-particle forces are added to the two- 
particle interaction in the total Hamiltonian. The exception is Eq. ()2.25|) for the mass 
operator in which the many-particle forces enter in the explicit form. This result of the 
GGFF could be expected, but it needed to be proved. It allows to extend the standard 
GF methods developed for the Fermi systems with two-particle interaction to the systems 
interacting through the many-particle effective forces. 
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In the final equations we can set £ = 0. So, in what follows we shall omit the ^-indices 
of the functions implying the limit at £ = and coming back to the GFs without source 
field defined by Eq. (|2.12j) . The independent functional variable is the single-particle GF G 
satisfying the Dyson equation which follows from Eq. ()2.23j) : 

G(z u z 2 ) = G°(z u z 2 ) + J dz 3 dz 4 G°(z u z 3 ) E(z 3 , z 4 ) G(z A , z 2 ) . (2.44) 

In case of the Fermi systems without pairing correlations, when the external anomalous pair 
potentials vanish in Eq. (|2.5p . the above equations can be reduced to ones for the normal 
components of GFs which satisfy the condition: 

k 

E(*n - *n) G • • • , z k ; z[, . . . , 4) = . (2.45) 

71=1 

These reduced equations will coincide with ones obtained in Ref. JE]- On the other hand, 
if we restrict consideration of the Fermi systems with pairing to the case of two-particle 
interaction in the Hamiltonian (|2.4jl . i.e. if we assume: 

H = H° + V {2 \ (2.46) 

the Eqs. lET!Z2|) . ({ZZ5) , dZHED , and (j2~4l^ will coincide up to the permutation of the argu- 
ments with the corresponding equations in Ref. |17j . 

2.3 Transformation of the basic equations 

It is easy to see that any amplitude X possessing the properties ()2.H7j) and (|2.H8jl can be 
represented in the form 

T(zi,z 2 ; z 3 ,z 4 ) = \ [U(z 1 ,z 2 ; z 3 , z 4 ) - U(z 2 ,Zi, z 3 ,z 4 )) , (2.47) 

where the amplitude U satisfies the equalities 

U{zx,z 2 ] z 3 ,z 4 ) =U(z 4 ,z 3 ; z 2 ,Zt) =U(z 2 ,z 1 ; z 4 ,z 3 ), (2.48) 

but generally U does not possess the property of antisymmetry (|2.37jl . So in the following 
we shall refer to the U as the non-antisymmetric amplitude of the effective interaction. 
Obviously, the Eq. ()2.47|) does not define U uniquely. To specify the definition let us note 
that Eqs. ()2.37|) are fulfilled if the functional variable of £ | in Eq. ()2.36|) is taken in the form 
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()2.33|) . However, one can formally define a functional derivative in which the substitution 
()2.33|) is not supposed. Following this definition we set: 

U( Zl ,z 2 ;z 3 ,z 4 ) = i S ^ Zl \ , (2.49) 
dG(z 4 ,z 3 ) 

where the condition ()2.14j) is not supposed to be fulfilled under variations of G. 

As an example, let us consider the first order in the two-particle interaction. From 
Eqs. (J2~23|) . (l2~22l . and (J2~2T| we have: 

S [11 («i,%)=i / dz 3 dz 4 [W (2) (z u z 4 ; z 3 ,z 2 )-\W {2 \z 1} z 2] z 3 ,z 4 )}G(z 3 ,z 4 ). (2.50) 



£ M is the generalized Hartree-Fock-Bogoliubov (HFB) contribution into the mass operator 
[the term generalized is used since the exact single-particle GF enters Eq. (|2.5()j) ]. Using, 
further, the definitions ()2.49|) and (j2.36J) we obtain from Eq. ()2.50|) : 

UW( Zl ,z 2 ; z 3 ,z 4 ) = W (2) (z 2 ,z 3 ; z u z 4 ) + lW (2 \z 2 ,z i; z 4 ,z 3 ), (2.51) 
l [1] ( Zl ,z 2 ;z 3 ,z 4 ) = \[uW{ Zl ,z 2 - z 3 ,z 4 )-UW(z 2 , Zl ; z 3 ,z 4 )}. (2.52) 

It can be seen that Eqs. (J2.48j) for the U W follow from the antisymmetry of the W ( k > and 
from the property (j2.30J) . while Eqs. (j2.37|) are not fulfilled. 

For the further applications it is convenient to introduce a non-antisymmetric response 
function R satisfying the following BSE: 

R(zx, z 2 ; z 3 , z 4 ) = R°( Zl ,z 2 ; z 3 ,z 4 ) + i J ' dz 5 dz 6 dz 7 dz$ 

x R°(z 1 ,z 2 ; z 5 ,z & )U(z b ,z & ] z 7 , z 8 ) R(z 7 , z s ; z 3 ,z 4 ), (2.53) 

where 

R°{ Zl , z 2 - z 3 , z 4 ) = -G(z 3 , Zl ) G(z 2 , z 4 ) . (2.54) 
It is easy to show, first that the following equalities hold: 

R{zi, z 2 ; z 3 , z 4 ) = R(z 4 , z 3 ; z 2 , zi) = R(z 2 , Zl ; z 4 , z 3 ) , (2.55) 

while the Eqs. f!2.39|) do not hold for R. Second, the response function L, which satisfies 
Eq. ()2.41j) . is expressed in terms of R as: 

L(z 1 ,z 2 ] z 3 , z 4 ) = R(zt, z 2 ; z 3 , z 4 ) - R(z 2 , Z\\ z 3 , z 4 ) . (2.56) 

Thus, for the determination of the response function L it is sufficient to solve Eq. ()2.53|) for 
the non-antisymmetric function R. 
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Now, following the method described in Ref. jSj, we represent the total mass operator 
£ and the total non-antisymmetric amplitude of the effective interaction U as a sum of two 
terms: 

£ = £ + £ e , U = U + U e , (2.57) 

where 

£(z u z 2 ) = S(y 1; y 2 ) S(ti — t 2 ) , (2.58) 
U(zi, z 2 ; z 3 ,z 4 ) = U(yi,y 2 ; V3,y^) S(t 1 -t 2 ) S(t 3 -t 4 ) S(ti-t 3 ) . (2.59) 

After transformation to the energy representation (see Sec. l2.4|) the first terms in Eqs. (|2.57|) . 
£ and U, are found to be energy- independent. The term £ corresponds to the mean-field 
contribution into the mass operator including the pair potentials. The term IA corresponds 
to the residual energy-independent interaction both in the ph, and in the pp channels. The 
second terms in Eqs. ()2.57|) . S e and U e 1 have a strong energy dependence and represent 
dynamic contributions of complex configurations. 

We stress that the quantities £ and U (and, consequently, £ e and U e ) are not defined 
rigorously by Eqs. ()2.57)) . They will be specified in the following within the framework of the 
model to be considered. At the moment only the general properties, which are expressed by 
Eqs. (j2.58|) and ()2.59|) . are important. Notice, however, that in particular case of the self- 
consistent HFB approximation restricted by the two-particle interaction we have: £ = E M , 
U = U^ l \ where the right-hand sides are defined by Eqs. ()2.5()j) and (|2.51j) with the exact 
GF G being replaced by the HFB GF G which is the solution of Eq. (|2.61|) (see below). 

Using decompositions ()2.57|) one can transform both of Eqs. ()2.44j) and ()2.53|) to the 
system of two equations. In the symbolic notations we have: 



G = G + GT l e G , (2.60) 

G = G° + G°tG, (2.61) 

R = R e + iR e UR, (2.62) 

R e = R° + iR°U e R e . (2.63) 



Proceeding by the same method as in Refs. [HI [TU] the last equation can be brought to the 
following form: 

R e (z 1 ,z 2 ; z 3 ,Z4j = R°(zi,z 2 ; z 3 ,z 4 ) + i J dz 5 dz 6 dz 7 dz 8 

x R°(zi, z 2 ; z 5 , z & ) W e (z 5 , z 6 ; z 7 , z 8 ) R e (z 7 , z 8 ; z 3 , z A ) , (2.64) 
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where 

R°(z 1 ,z 2 ; z 3 ,z 4 ) = -G(z 3 , z x ) G(z 2 , z 4 ) , (2.65) 
W e {z h z 2 ; z 3 ,z 4 ) = V e (z 1 ,z 2 ; z 3 , z 4 ) - iZ e (z 3 , z 1 )Z e (z 2 , z 4 ) , (2.66) 

, Z 2 \ Z 3 , Z4) 1 z 2 ] Z 3 , Z4) 

+ iS e (2 3 , zi) G~ 1 (z 2 , z 4 ) + iG~ l (z 3 ,z 1 ) T l e (z 2 , z 4 ) . (2.67) 

The Eq. (f2~64jl is the basic one for the building of model which will be considered in the 
second part of the paper. 

2.4 Single-quasiparticle basis functions and the energy representation 

For the following analysis it is required to introduce a set of basis functions {V ; 1 (y)} in the 
extended space defined previously in Sec. 12.11 The usual conditions of orthonormality and 
completeness are supposed to be fulfilled: 

I dyPMMv) = Kv , E ri(v)A(v') = <W) • ( 2 -68) 
J 1 

It is convenient to consider ip^y) to be the eigenf unctions of the operator: 

H(y,y')=H°(y,y') + t(y,y'), (2.69) 

where 7i° defines the single-particle term of the total Hamiltonian according to Eqs. ()2.9|) 
and (|2.1()jl . S is the mean-field contribution into the total mass operator in Eqs. ()2.57jl . 
Thus, we assume the following equation to be fulfilled: 



dy"H(y,y , )4> 1 (y') = E 1 4> 1 (y). (2.70) 

Since 7i possesses the same symmetry properties as the operators 7i° and E, i.e.: 

H(y,y / )=H*(y',y) = -H(y',y), (2.71) 

it is not difficult to see that the complete set of the eigenfunctions of 7i is divided into two 
equal parts which are related by the operation of conjugation: 

My) = rM- (2.72) 

For the corresponding eigenvalues we have: E\ = —Ex- So one can denote: 1 = {Ai.,771}, 
1 = {Ai,— T]i}, where Ai is the index of the usual single-particle configuration space (e.g., 
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A = {t\, n, I, j, m} for the spherically symmetric system), rji = ±1 is the sign of the eigenvalue 
E r . 

E 1 = ij 1 E Xl , £ Al = |#il- (2.73) 
In the representation of ^-functions, the operators b defined by Eq. (j2.1j) have the 

form 

6 1= f dyM{y)b(y). (2.74) 



It is worth noting that the 6-operators in this representation are simply related to the creation 
and annihilation operators of the quasiparticles a\ and a x which are usually introduced in 
the HFB theory. Namely, we have (see Ref. [Ej): b x+ = a x , b x _ = a x . So, in what follows 
we shall refer to the functions ip^y) as the single-quasiparticle functions. 

A more detailed form of the functions ip-^y) can be obtained making use of the Bloch- 
Messiah theorem, see Refs. To formulate the result let us note that for time- reversal 

invariant fermion system, the set of the single-particle indices A can be divided into two 
subsets of conjugate indices p and p which represent "paired" states, i.e.: {A} = {p} U {p}. 
In particular, for spherically symmetric system the conjugate indices are A = {(A),m} and 
A = {(A), — m} where (A) = {t\, n, l,j}. According to the Bloch-Messiah theorem the single- 
quasiparticle functions can be represented in the following form (see Ref. 



A' A' 

where C x , x is a unitary matrix and the functions tp\ )V (y) have the form: 



4> Pt +(x,-) = -v p (p*(x), tp P)+ (x,-) = v p <f* p (x), 

Ip p ,-(X,+) = -V p lfi p (x), Tp p -(X,+) = V p if p (x), 

ip p -(x,-) = u p ip* p (x), tp P: -(x } -) = U p if* p {x). 



(2.76) 



Here {ip\(x)} is a complete set of orthonormal functions in the usual single-particle space, v\ 
and u\ are real non-negative numbers which satisfy the following conditions: u x = y/l — v x , 
v p = Vp ^ 1. For the spherically symmetric system we have: 



[2.77) 



$a,+(z,+) =u x ip x (x) ) ^a,+ (x,-) = (-l) l+j+m v x v* x (x), 
4>x,-(x,-) =u x <p* x {x), $x,-(x,+) = (-l) l+ i +m v x <p- x (x). 

Up to now we did not restrict our analysis to the systems where the number of 
particles is conserved exactly. The reason is that our aim was modification of the existing 
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general GF formalism for the arbitrary Fermi systems with pairing. However, application of 
the formalism to the atomic nuclei, we are interested in, implies that the particle-number 
conservation law is fulfilled. Thus, in what follows we assume that the total Hamiltonian 
H defined by Eqs. (|2.4j) - (|2.7|) does not contain the external anomalous pair potentials A 
and that Eqs. (|2.8J) hold. In that case the ground-state wave function | ) which enters 
the definition of the GFs (j2.12|) is an eigenfunction of the particle-number operator. This 
means that the exact GFs do not contain anomalous components. In particular, the exact 
single-particle GF satisfies condition [cf. Eq. (|2.45j) ]: 

( X -x')G(z,z') = 0. (2.78) 

However, this condition is not fulfilled for the GF G which is the solution of Eq. (|2.61|) 
with the operator E including the pair potentials independently of the Hamiltonian H (e.g., 
within the HFB approximation). Justification of using such GF G is as follows. It enables 
one to take into account pairing correlations effectively and should be considered only as 
an approximation to the exact GF. The latter is found from Eq. (12.60(1 in which the mass 
operator S e has to contain all necessary corrections to S, such that the solution of Eq. (|2.fi()j) 
satisfies Eq. (12.78(1 . Of course this scheme should be considered only as a philosophy of the 
approach, i.e. as an ideal program which is difficult to implement completely in practice. 

Let us now define the energy representation of the Green functions and of the related 
quantities entering above equations. Making use of the basis (y)}, let us introduce the 
following Fourier transformations: 

G?ia(e) = dz 1 dz 2 il>l(y 1 )i/; 2 (y2)5(t 2 ) ex.p(ie(t 1 -t 2 ))G(z 1 ,z 2 ), (2.79) 



R\2,za{ u ) = -i \ dzxdz 2 dz 3 dz^ipy (yi) ip 2 (y 2 ) ip 3 (y 3 ) V>4 (2/4) 

x 6(ti -h-0) 5(U - t 3 - 0) 5(U) exp(iuj (t 3 - h) ) 

x R(z x ,z 2 ; z 3 ,z 4 ) , (2.80) 

Wi 2 ,34(<*>, £, e') = J dz x dz 2 dz 3 dz 4 Vi(yi) ^2(2/2) ^3(2/3) ^{Ui} 

x 6(U) exp(icu (t 3 - ti) + ie (t 2 - ti) + ie' (t 3 - t A ) ) 

X U{z!,z 2 ; z 3 ,z 4 ) . (2.81) 

The quantities £r 12(e), E 12 (e), Ef 2 (e), ^1234^)) ^1234(^5 e i e ')i an d others are defined in 
analogy to these formulas. In accordance with Eqs. (|2.56| . ()2.72|1 . and (j2.80() we have: 

L 12M {uj) = #12,34(0;) - #21,34^) • (2.82) 
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Notice that spectral expansion for the response function L(u>) has the form: 

W) = "EE r P " ■ m , (2-83) 

which is similar to the analogous formula for the Fermi systems without pairing correlations 
(see, e.g., Ref [IS!)- The difference consists in the definition of the transition amplitudes. In 
Eq. QZBfl we have: 

P n i2 ] = K+i (n\b\b 2 \0) + 5 Vj _ x (0\b\b 2 \n), oo n = E n - E , (2.84) 

where | n ), | ), E n , and E Q are the eigenf unctions and the eigenvalues of the Hamiltonian H . 
Taking into account the definitions (|2.1j) and (|2.74jl . one can see that even if the ground-state 
wave function | ) is an eigenfunction of the particle-number operator, the amplitudes p™ 2 
take nonzero values not only for the transitions between the states with the same number of 
particles, but also for the transitions between the ground state of the iV-particle system | ) 
and the states | n ) of the systems consisting of iV ± 2 particles. Thus, the spectral expansion 
(12.83)1 contains information about excitations of the iV-particle system both in the ph, and 
in the pp channels. 

3 QUASIPARTICLE TIME BLOCKING APPROXIMATION 
3.1 General framework 

Let us now turn to the question of determining the physical observables and related quantities 
in this approach, namely, excitation energies u n and transition amplitudes p™ 2 ■ It follows 
from the Eqs. ()2.82|) and (j2.83J) that in order to find these characteristics we need to know 
the response function R(u>), i.e. to solve the system of equations ()2.62|) and ()2.64|) . The basic 
difficulty in solving this task is that the Eq. ()2.fi4j) contains energy- dependent interaction U e 
and mass operator S e . Notice, however, that these energy-dependent quantities arise only 
in the case when the dynamic contributions of complex configurations are taken into account 
explicitly. In the energy representation Eq. ()2.64|) is an integral equation for the function 
R e (u, e, e') [defined in analogy to Eq. ()2.81|l ] over the energy variable e, which cannot be, 
strictly speaking, reduced to the closed equation for R e {ui) due to the energy dependence of 
U e and S e . Fortunately, there are methods that allow us to avoid complicated problem of 
exact solution of this equation, making use of certain approximations. One of such methods 
will be considered here. 
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We begin by noting that if we use the eigenfunctions of the operator H, i.e. the set 
{ifrAy)}, as the basis functions, the single-particle GF G is diagonal: 

G 12 (e) = J" . , (3.1) 

e - E x + z -0 

as follows from Eqs. (|2.24jl . (|2.filj) . (|2.fi9|) . and (|2.7()jl . In the time representation we have: 

Gi 2 (h, t 2 ) = -i r]i5 12 0(7] 1 t 12 ) exp(-iEiti 2 ) , (3.2) 

where t\ 2 = t± — t 2 , 9(r) is the step function. These expressions are formally identical 
with analogous formulas for the normal GFs, except that they are written in the extended 
basis representation. It enables one to apply, practically without changes, the method of 
chronological decoupling of diagrams (MCDD) to the solution of Eqs. ()2.62|) and ()2.64|) 
which contain the GFs with pairing. The MCDD was developed in Ref. ^0j for the solution 
of the ph-channel BSE in the normal Fermi system including dynamic effects both in the 
interaction and in the mass operator. The idea of the method is similar to that used in 
the other methods developed earlier for the solution of the analogous problems in Refs. [201 
(ph-channel BSE) and [2U E2J ES] (pp-channel BSE). However, the MCDD differs from the 
aforementioned methods in some details, in particular concerning the treatment of the GSC2. 
Almost all the resulting equations obtained by means of this straightforward extension of 
the MCDD are found to be formally identical with equations for the normal Fermi system 
in the same sense as Eqs. (13 .1|) and ()3.2j) . Because derivation of these equations in the latter 
case was described in detail in Refs. [31 EH], we will draw only the main formulas and the 
final results. 

First of all, the function R° — —GG entering Eq. ()2.64|) is divided into two parts: 
£o = £ o(a) + ft o(b) where 

^12,34 (*i> t 2 ; £3,^4) = —5 Vi _ V2 6(i]it 41 ) 6(r]it 32 ) G 3 i(t3,ii) G 2 4(t 2 ,U) , (3.3) 

R°( b ' is the remainder term, which is absorbed in part in the renormalization procedure. 
As compared with the initial function R° defined by Eq. (|2.65|) . the term R ^ contains 
two additional time-dependent step functions and the factor S v _ v which play twofold role. 
On the one hand, they allow to obtain closed set of the algebraic equations in the energy 
representation for the main component of the function R e {ui) [see Eqs. ()3.5|) . (j3.17p . ()3.18|) 
below] which is much more simple for the solution as compared with the initial Eq. ()2.64[) . 
On the other hand, owing to these additional 9 and 5 functions, an approximate solution 
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of the Eq. ()2.64j) obtained by this way contains the contributions of the 2q and 2qCS>phonon 
configurations, while more complicated intermediate states (e.g., 2q®2phonon, 2qCg>3phonon, 
and so on) are blocked, in part, in the time representation. So in the following this scheme 
will be referred to as the quasiparticle time blocking approximation (QTBA). 

Further, renormalization procedure is applied to Eq. (j2.62j) for the response function 
R, which leads to the following equation for the effective response function R in the energy 
representation: 

^34 M = Au M (v) ~ ^12,56^) -^56,78 ^34 M , (3-4) 

5678 

where A(u) is a joint (ph and pp) correlated propagator, T is an amplitude of the effective 
interaction. The propagator A{lo) is the main term of the formal decomposition 

R$2 M (w) = A 12M {w) + S 12i34 . (3.5) 

It contains: (i) the sum of an infinite number of terms to all orders in R°( a \ and (ii) the terms 
linear and quadratic in R ^ which are related to the GSC2. The term B in Eq. ()3.5|) is an 
auxiliary quantity which is supposed to be energy-independent. In addition, it is supposed 
that B contains all the contributions which are not included explicitly in the propagator 
A(u). The effective interaction T and the effective charge operator e are defined by the 
formulas: 

•^12,34 = e 12, 56 ^56,34 ' e 12,34 = ^13 ^24 ~ •^12,56 -^56,34 = ( e 34,12) • (3-6) 

56 56 

In terms of these quantities the exact response function R is related to the effective response 
function by the ansatz: 

-^12,34 ( u ) = e 12,56 -^5^78 C^) e 78,34 + ^12,56 e 56,34 ■ (3-7) 
5678 56 

One of the basic quantities, which determines the physical observables in this ap- 
proach, is the nuclear polarizability II(a;). More precisely, it determines distribution of the 
transition strength caused by an external field V°(x, x'). The function II (u;) is defined as 

nH = -\ £( e O*i<34M (ey%, ( e y°) 12 = £ e 2M3 v& , (3.8) 

1234 34 

where 

V° = J dydy'^M^y')^ [6 Xl+ iV°(x,x') - 5 x ^V°(x',x)] = -V£. (3.9) 
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In particular, the strength function S(E) which is frequently used for the description of 
nuclear excitations is expressed in terms of the polarizability as 

S(E) = -L Im V V 2 °* R 12M (E + i A) = -- Im IL(E + i A) , (3.10) 

where A is a smearing parameter. The formulas ()3.4j) - (j3.8j) . ()3.1()|) are completely analogous 
to the ones for the normal Fermi system (see Ref. [Sj), except for the factors \ in Eqs. (|3.8j) 
and ()3.10|) which arise due to definition ()3.9|) of the operator V in the extended space taken 
in the antisymmetric form. 

3.2 Correlated propagator within the QTBA 

Eq. (|3.4j) for the response function is still quite general. To formulate a model we have to 
define the correlated propagator A(uS). In particular, if we neglect the dynamic contributions 
of complex configurations, i.e. if we put U e = and S e = 0, we come to the QRPA. In this 
case we have: A(u) = A(u) where A(u) is uncorrected QRPA propagator: 

4^(") = - ^'^i U * M . E a = E 1 -E 2 . (3.11) 

To go beyond the QRPA we have to find reasonable approximations for the quantities lA e 
and S e . In the present work we shall use a QRPA-based version of the quasiparticle-phonon 
coupling (QPC) model (see Ref. [24"]). This model is discussed and used in a variety of 
papers: see, e.g., Refs. j2 El HI EE3 UH| ■ Within the QPC model one can restrict oneself to 
so-called g 2 approximation where g is an amplitude of the quasiparticle-phonon interaction 
(see Ref. jS] for more details). Under some simplifying assumptions, this approximation 
can be obtained in the GF method (see Refs. |23 EE1 Elj ) • Within the QPC model and g 2 
approximation we have the following formulas for the quantities U e and S e : 



r],m 



e - e' + rj (u m - i-0) ' 

e 771(77)* m(ri) 

°V,V 3 9 13 923 



12<) &'-*->>(<*.-*■<»• ' ' 

where 77 = ±1 . Hereafter it is assumed that the quasiparticle-phonon amplitudes <fo are 
related to the transition amplitudes Pi 2 by means of QRPA equations: 

m(rj) _ \ ^ r m(rj) m(rf) _ T h Vir-Jh m(v) (o -ia\ 

9 12 — 7, -M2.34 P34 > Pl2 ~~ P~ 9l2 : l* 3 - 1 ^ 

^ V^m — &12 
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where T is an amplitude of the effective interaction which generally differs from the amplitude 
T entering Eq. (|3.4jl . Notice that the QRPA equations acquire very simple form in the 
representation of single-quasiparticle -0-functions. The Eqs. ()3.14j) have to be supplemented 
by the normalization condition 

^£^iK 2 ( T = i (3.i5) 

12 

and by the condition of the antisymmetry 

= _ p ^ V ) (3 16) 

which is obviously fulfilled for the exact transition amplitudes defined by Eq. ()2.84|) . Notice, 
however, that in contrast to the case considered in Refs. [3J EI] , the quasiparticle-phonon 
amplitudes in Eqs. (13.12)1 and (|3.13|) determine the coupling with excitations both in 
the ph, and in the pp channels [see Eqs. ()3.14j) and comments after Eq. (|2.84j) ]. 

The QTBA in combination with Eqs. ()3.12|) and f)3. 13j) leads to the following ansatz 
for the correlated propagator: 

A2, 34 M = E [*w 6 ™ + Q££(P>) ] 4e,^H [5 73 5 8i + ] + P£#(w) , (3.17) 

5678 

where the upper indices denote products of the first and the second pairs of lower //-indices. 
In particular, for the component A^^^u) of the propagator we have: = r] 3 r] 4 — — 1. 

This component is determined by the equation 

Au&iw) = ^12,34 M -E^H $56,78 M 4^34 M , (3.18) 

5678 

where A{u) is the QRPA propagator defined by Eq. (|3.11)1 . For the remaining quantities in 
the Eqs. (l3~T7j) and (l3~TH|) we obtain 

QiSiH = Qitd^ + ^J^J ^S 24 - 5 31 ^) , (3.19) 



VGSC 
^31 


^24 


- 5 31 


VGSC 
^ 24 




-^24 




VGSC 
^31 


^24 


- <*31 


VGSC 
^24 


^31 


-^24 





Q&JiH = Q { uPr^)-K,-r,J V3>Vi [^S u -5 31 ^) , (3.20) 

$12,34^) = $lt34M + $i?34 + $l G 2 S 34 S - e -M- (3-21) 

In these formulas, a superscript "res" denotes the resonant parts of the amplitudes, the 
quantities $ GSC and $ GSCs - e (u;) represent contributions of the GSC They consist of the 
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static part arising from the induced interaction ($ GSC ) and of the part arising from the 
self-energy insertions ($ GSC s e -): 

*S£ = ^V^^r^Ekfe^f^^fc^fl, (3-22) 

77, m ^ ' 

*i G 2 S 34 S - e -H = ^i^,-, 2 ^3,-, 4 (sS sc (5 24 + g 24 )-(5 3 i + 93i)S^ 4 sc 

- ( 931 $24 + $31 924 + 931 9 24 ) [ <*> ~ \ ( E 12 + E34)] J , (3.23) 

where 

£rf c = + ^^^J ^ 77 ^ )7 , 3 (pS (T?) *^ (T7) + ^ (77) *^ (r?) ) , (3-24) 

3,r/,m ^ ' 

9i 2 = E <W • (3-25) 

3,r],m 

The component P( ++ \u) of the correlated propagator and the resonant parts of the ampli- 
tudes entering Eqs. (|3.19j) - (j3.21j) are defined as 



p(++) 

-"12,34 1 


» = 


V^ 1 >m56(?j) .(r;) 
/ y S12 ^56,78 1 
5678, 77,771 


[uj 




> 777.78(77)* 
S34 > 


(3.26) 


n (+-)res 
Vl2,34 1 


M = 


A"i56(»y) 7(77) 
/ y ^12 ^56,78 1 
5678,77,m 




-r)u m ) 


m78(r ? )* 
734 J 


(3.27) 


n (-+)res 
Vl2,34 ' 


» = 


/ j /12 ^56,78 1 
5678, 77,771 






,m78(7?)* 
1 S34 ) 


(3.28) 


fates 
^12,34' 


» = 


ru m56 ( r ?) 2(»7) , 
/ y '12 ^56,78 1 




- T)UJ m ) 


77178(77)* 
1 734 ) 


(3.29) 



5678,r;,m 



where 



= $ v , % K,-v 2 $v 5 ,-v 6 ( *i5 9$ V) ~ 9? 5 { %2 ) , (3-30) 

Cl2 ^ = ^»?,r; 5 ^,7/2 $Vs>-V e ( ^ 62?7 _ Pl5^$62 ) 5 (3.31) 

A( + )(u;) and AH^) are the positive and the negative frequency parts of the QRPA propa- 
gator defined by Eq. flUTE), i.e.: A(u) = i (+) (tu) + > 

Correlated propagator A (a;) defined by Eq. (|3.17|) and subsequent equations includes 
contributions of three types: (i) pure 2q configurations associated with uncorrelated QRPA 
propagator A(u), (ii) 2qCg>phonon configurations introduced by the quantities U e and S e , 
and (iii) uncontrollable more complicated configurations arising due to the GSC effects and 
their combinations with the above-mentioned configurations. 
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3.3 Sum rule analysis and a refinement of the model 



The formulas of the previous subsection completely determine the correlated propagator of 
the model within the g 2 approximation. By construction, this propagator contains all the 
g 2 contributions, including those from GSC However, exact fulfillment of the sum rules in 
this approach is not guaranteed. Let us consider this question in more detail. Usually, the 
sum rule is meant in the sense of relation between the moment m k of the strength function 
S(E) and the ground state expectation value of certain operator (see, e.g., Ref. |2B])- The 
moment is defined as 

1 f°° 

m k = - S(E) E k dE (3.33) 

^ J-oo 

at A — > +0 in Eq. (|3.1Uj) . Introducing asymptotic expansion of the exact response function: 



(fc) ,-fc-l 

1V1 

and using Eqs. (|2.82|) . (|2.83|) . (|3.9|) . (|3.1U|) . one can show that the moments m k are expressed 

(k) 

through the coefficient functions M 12 34 by the formula: 

mfc = lE y 2 i*M 1 g 4 V r 4 3 . (3.35) 

1234 

In particular, making use of the BSE ()2.53j) in the energy representation one can obtain: 

= <*31 P24 - PSI ^24 , (3-36) 

where p 12 = ( | b\ b 1 | ) is the extended density matrix (EDM). Substituting Eq. ()3.36|) 
into the Eq. ()3.35|) we get the so-called non-energy- weighted sum rule (NEWSR): 

m = lTr(p[V°,V^}). (3.37) 

Notice that the factor j in this formula (instead of the usual factor |) arises from the 
definition (j3.9|) of the external field operator in the extended space taken in the antisymmetric 
form. 

Thus, in order to ensure exact fulfillment of the NEWSR the coefficient function 
M/2 34 of the model must have the form ()3.36|) with properly normalized EDM p. Using the 
formulas of the previous subsection and the definition of the quantity 34 through the 
expansion ()3.34|) one can show that Eq. ()3.36|) is not fulfilled in the model considered and 
that the NEWSR is fulfilled only up to within the terms of order g A . 
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It is not difficult, however, to remedy this drawback within the above- described 
scheme based on the g 2 approximation. First of all, let us include energy-independent op- 
erator E GSC defined by Eq. (J3.24)) into the mean-field part S of the total mass operator 
It can be done because the only constraint on the operator £ was the condition of its 
energy independence. This redefinition of S means that we have to use in all the equations 
the operator S e (e) = ^ e {s) — £ GSC instead of H e (e). The replacement of £ e (e) by S e (e) 
leads to disappearance of all the terms containing £ GSC in Eqs. (J3.19)) . 1)3.20)1 . and ()3.23j) . 
The remaining part of the amplitude $ GSC s - e, {uj) in Eq. (j3.21)l can be taken into account 
through the renormalization of the QRPA propagator A(u) within the g 2 approximation. 

To this aim let us introduce matrix ^12,34 defined by the following equations: 

2^ ^12,56 "^56,34 = ^31 ^24 — <?31 ^24 — ^31 ?24 ; ^12,34 = ^34,12 ) (3.38) 

56 

where the matrix q 12 is defined by Eq. ()3.25j) . In addition it will be supposed that the matrix 
^12,34 is positive-definite that can be always fulfilled if all the eigenvalues of the matrix q 12 
satisfy condition: q%<\ [notice that ^ 0, as follows from Eq. ()3.25j) ]. Because according 
to Eqs. (j3.14j) and (J3.25)) we have: q 12 = 0(g 2 ) , the pointed condition is consistent with the 
previous model assumptions. Thus, from Eq. (J3.38)) it follows that 

^12,34 = 4l S 24 ~ | ( ?31 ^24 + <*31 ?24 ) + 0(g A ) . (3.39) 

Further, using Eqs. (j3.ll)) . ()3.23j) (without terms containing operator S GSC ), and 
()3.39)1 . we obtain: 

A(u) - A(uo) $ GSC s - e \u) A(uo) = Z A(u) Z + 0(g 4 ) . (3.40) 
It enables one to redefine correlated propagator replacing Eq. (j3.17j) by the following ansatz: 

AumM = E Z k™^) 4 6 ,78^) + P£${u) , (3.41) 

5678 

where 

^2,34^) = E^IS^ + ^SS^M] ^56,34, (3-42) 
56 

^2,34^) = E^,56^53^64 + Q56,3i reS M]- (3-43) 
56 
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In Eq. (13.41)1 the propagator (u) is determined by Eq. (|3.18j) in which the amplitude 

$(o;) is now defined as 

®12,sM = E ^12,56 [ + ^5 G 6 S 78 ] ^78,34 , (3-44) 

5678 

instead of the Eq. (jS2H)- 

It is easy to see that the propagator A(u) defined by Eqs. ([3~T8]l . ([3~3S]l . (jS3H)-(|S31|) 

coincides up to within the terms of order g A with the propagator defined in the previous 
subsection [see Eq. ([3.17)1 and subsequent equations]. On the other hand, assuming that 
the effective charge in Eq. ()3.7|) is equal to unit operator and making use of the expansion 
([3.34)1 . one can find that this new redefined propagator leads to the following result for the 
coefficient function M^u- 

M 12°34 = Vl V - V2 ( *31 6 24 ~ 931 5 24 ~ <*31 ?24 ) + Vl \ , r? 2 ( 5 31 ?24 ~ 9 3 1 5 24 ) ■ (3-45) 

Here the first term containing $ v t - v follows from Eq. (|3.38j) . The second term arises from 
the contribution of the component P( ++ )(u;) in Eq. (|3.41)1 . From the Eq. (J3.45)l we obtain 
that in the modified version of the model the coefficient function M/2 34 has the form (J3.36)l 
with correlated EDM p defined as 

P12 = P12 + Vi 9i2 > (3-46) 

where p 12 = S v _ 1 8 12 is the EDM of the HFB theory in the representation of ^-functions. 
Therefore we conclude that if the EDM ([3.46)1 is normalized by the usual condition 

I dy>5 x ^ +l 5 T ,^^i(y f )r2(y')Pi2 = N T , (3.47) 
J 12 

the NEWSR is fulfilled exactly within the QTBA. It is worthwhile mentioning that the EDM 
([3.46)1 arising in the QTBA coincides with the correlated EDM p c which can be obtained 
from the Dyson equation ([2.60)1 and the Eqs. ([3.1)1 . ([3.13)1 . ([3.24)1 within g 2 approximation: 

f°° dr 

P 12 = Pi2 = J ^e^[G(6) + G(e){^(e)-Z GSC )G(e)} 12 , r^+0. (3.48) 

Finally, note that the GSC effects included by means of renormalization of the QRPA 
propagator within the QTBA with the help of the matrix Z are the same as the effects 
included in the renormalized QRPA (RQRPA, see, e.g., Ref. 6 ). It is known (see Ref. [25] ) 
that within the standard RQRPA, the Ikeda sum rule (being a particular case of the NEWSR) 
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is violated. The above analysis allows us to understand the reason of this violation. It follows 
from Eq. ()3.45|) that to satisfy the NEWSR within the QTBA it is necessary to take into 
account contribution of the term p( ++ \u). This term represents dynamic contributions of 
the GSC into the QTBA propagator (|3.41J) . which cannot be reduced to the renormalization 
of the QRPA propagator A{oS) and which are absent in the RQRPA. 

3.4 Antisymmetrization of the equations 

and inclusion of the two-phonon configurations 

As can be seen from the spectral expansion ()2.83j) . the physical observables of the theory 
are completely determined by the antisymmetric response function L(u). Within the QTBA 
this exact function is approximated by the effective antisymmetric response function L (uj) 
defined as [cf. Eq. fl2~%2jl ] 

Lgs^w) = RSjuifj) - J2g M (uO , (3.49) 

where R eS (uj) is the solution of Eq. (|3.4j) . It is easy to prove the following. First, the polariz- 
ability n(u;) [see Eqs. (|3.8j) ] is actually determined by the function L (u), while symmetric 
part of the function R eS (uj) does not contribute to the Eqs. ()3.8|) . In other words, II (u;) is 
invariant under the transformation: R eS (uj) —>■ ^L cS (uj). Second, the function L eG (u) is the 
doubled solution of the antisymmetrized QTBA equation obtained from the Eq. ()3.4j) with 
the help of antisymmetrization of the correlated propagator A(u). This antisymmetrization 
can be implemented by means of the following transformations in Eqs. (j3.18|) . ()3.26|) - ()3.29|) : 
A{uj) — > ^L°(uj), A^'ioS) — > | L ^ (to), where L°(uj) is the antisymmetric (uncorrelated) 
QRPA propagator, L ot ~ :r, \oj) represents its positive and negative frequency parts [cf. Eqs. 
(ETTTJ) and (EH2J) ]: 

VLmM = £ , = (^3 «^ - ^3 . (3.50) 

We did not use the antisymmetric form of the QTBA equations from the very beginning to 
simplify their derivation and analysis. However, the antisymmetrization facilitates numerical 
solution due to reduction of the dimensions of matrices entering these equations. 

The model described above allows for the following straightforward extension related 
to the definition of the resonant parts of the amplitudes entering Eqs. (|3.18jl . (|3.41jl — (|3.44|l for 
the correlated propagator of the QTBA. Contributions from these resonant parts [defined by 
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Eqs. fl3.26|) - ()3.29|) ] to the response function describe simultaneous propagation of the phonon 
and of the uncorrelated quasiparticle pair. Natural generalization of this model is inclusion of 
the correlations in the quasiparticle pair entering 2q®phonon configuration, i.e. replacement 
of the uncorrelated pair by the phonon. For the ph-channel BSE in the normal Fermi system 
similar generalization, corresponding to the replacement of the lplhCg>phonon configurations 
by the two-phonon intermediate states, was discussed in Ref. [3U]. For the pp-channel BSE 
analogous procedure has been implemented in Ref. |2*2*] . 

Within the QTBA, two-quasiparticle correlations in the 2q®phonon intermediate 
states (i.e., two-phonon configurations) can be included in the following way. Correlated 
counterpart of the above-defined quantity L°™(o;) is L™(u;) representing the positive and 
the negative frequency parts of the antisymmetric QRPA response function L(to): 

n{ri) n(?7)* 

= E > = - £ VPl2 PU • (3-51) 

— , , tO — 77 tO n 

n=±i n 1 n 

In the Eq. (|3.51|) it is supposed that the QRPA energies uo n and transition amplitudes p™2 
satisfy Eqs. ()3.14jl — (jH.lfij) . The above considerations imply that transition to the two- 
phonon configurations within the QTBA can be accomplished by means of the replacement 
A^\uj) -> \ L^\u) in Eqs. (j3~2l)|) - (l3~2l?|) . which leads to the following result: 

i >mn(tj) >m»(i))* 

(++)/, ,\ 1 ST- ^^12 C 



P 12 o 4 [LO 



1 „ /" "iii-VI) /- unity i) 

- y ^ ™ , (3.52) 

2 ^ to-rjto n 



v ,m,n ■ mn 

. mn(rj) mn(7j)^ 



Q[td ies n = -\ E 12 734 — > ( 3 - 53 ) 

2 *— ' tO — T]LO n 



. ■ / mn 

r),m,n 

mn(rj) j-mn(rj) 



= 4E " 712 34 > ( 3 - 54 ) 

2 ^— ' to — r\u3 n 



r),m,n 

mn(rj) mn{rj)^ 



i E ' 7l2 734 , (3.55) 



2 z — ' to — rito 

r],m,n 



where = w m + to n 



7 12 =Z^ 7 12 P56 > Cl2 U -2^Cl2 P56 ■ ( 3 -56) 

56 56 

Physical arguments in favor of using Eqs. ()3.52|) - (j3.55|) instead of ()3.26|) - ()3.29|) are 
clear. Notice, however, that the derivation of the Eqs. ()3.52j) - (j3.55|) has not been rigorous. 
It enables one only to assert that these formulas recover the original Eqs. (|3.26j) - (|3.29|) in 
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the limit of vanishing quasiparticle interaction. One can also show, using the completeness 
of the set of QRPA transition amplitudes, that the NEWSR is fulfilled exactly within two- 
phonon version of the QTBA. The rigorous derivation of the Eqs. (|3.52j) - (|3.55j) is based 
on the inclusion of the additional (third order in the quasiparticle interaction) contributions 
into the dynamic amplitude U e defined by Eq. (|3.12j) and will not be considered here. It is 
worth noting that inclusion of the two-phonon configurations in Eqs. ()3.52j) - (j3.55j) brings 
the model closer to the QPM Comparing the QTBA and the QPM one can infer that 
treatment of the GSC within the QTBA is more consistent. A more detailed comparison of 
these models is beyond the scope of the present paper. 

3.5 Self-consistent scheme 

Finally, we briefly outline the scheme which enables one to eliminate spurious states within 
the QTBA. These states, being a common problem of the microscopic theories, are associ- 
ated with existence of the non-trivial external field operators V° satisfying the condition: 
[H, V°] = [H°, V }, where H and H° are the total and the single-particle Hamiltonian, cor- 
respondingly [see Eq. (|2.4|) ]. Elimination of the spurious states within a consistent theory 
implies that they must have zero excitation energy. In terms of the GF method this means 
that the exact response function R(u) must have a pole at uj = corresponding to the spuri- 
ous states. In particular, it is well known that the QRPA response function R(u) satisfying 
the equation 

R(u) = A(u) - A{uj) F R{uj) (3.57) 

has such a pole at u) — if the interaction amplitude T is related to the mean-field operator 
X by the self-consistency condition 

T = btjb~p. (3.58) 

In the QTBA the situation is more complicated since the correlated propagator A(u) 
in Eq. (|3.4|) has no simple structure of the QRPA propagator A(u). To avoid this difficulty 
let us note that the exact solution of the Eq. ()3.4|) with the propagator A(cu) defined by Eqs. 
fl3.41|) . (j3.18j) can be represented in the form: 

R cS (uj) = [l-P( ++ \uo)F p (uj)] Z L {u)R eS {u)Z R {u) [I -F p (lu)P {++) (uj)] 

+ P (++) H - P (++) (cu) T p {u) P (++) (^) , (3.59) 
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where energy-dependent interaction amplitude JF p (c<j) and the renormalized response func- 
tion R eS (uu) satisfy the equations: 

jrP^ = f -FP(++)(uj)F p {uj), (3.60) 

R cS {u) = A(w)-i(u)f(w)^(w), (3.61) 

with 

F{u) = Z r (uj) F p (u) Z l (uj) + . (3.62) 

If the amplitude J-(0) coincides with the interaction amplitude T satisfying Eq. (j3.58j) . 
the renormalized response function R eS (u) has the pole at u = corresponding to the 
spurious states by the same reasons as the QRPA response function R(u). To ensure the 
fulfillment of the relationship .F(O) = T we use the fact that the interaction amplitude T 
entering Eq. ()3.4)1 has not been constrained so far by any conditions besides the property 
of its energy independence. Let us now assume that the amplitude T satisfies the following 
equation: 

T = F p + F p p( ++ \$)F , (3.63) 

where 

T p = [Z*(0)] -1 [.F-$(0)] [Z L {0)}~\ (3.64) 

and the amplitude T is determined by Eq. ()3.58|) . If the Eqs. ()3.63|) and ()3.64|) are fulfilled, 
then it follows from Eqs. ()3.6()|) and ()3.62j) that T($$) = T . Consequently both function 
R eS (u) and function R eS (uj) have the poles at u = corresponding to the spurious states. It 
means that these states are eliminated, at least energetically, within self-consistent version 
of the QTBA defined by the above equations. 

4 CONCLUSIONS 

In this paper the problem of the microscopic description of excited states of the even-even 
open-shell atomic nuclei is considered. The generalized Green function formalism (GGFF) 
has been presented and used to formulate the model including pairing, two-quasiparticle 
(2q), and more complex, first of all quasiparticle-phonon correlations. The GGFF is a 
modification of the existing versions of Green function formalism which is more suitable 
for solving the problem considered here. Within the GGFF the normal and the anomalous 
Green functions in the Fermi systems with pairing are treated in a unified way in terms 
of the components of generalized Green functions in a doubled space. This treatment is 
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analogous to the method used in Ref. ^7j. In the GGFF this method is extended to the Fermi 
systems interacting through the two-, three-, and other many-particle effective forces, that is 
of importance for the nuclear physics where the many-particle forces play an essential role. 
Within the framework of this formalism the generalization of the model of Ref. including 
the pairing correlations has been developed. The physical content of the model is determined 
by the quasiparticle time blocking approximation (QTBA) which allows one to keep the 
contributions of the 2q and 2q®phonon configurations, while excluding (blocking) more 
complicated intermediate states. It has been shown that within the QTBA the non-energy- 
weighted sum rule is fulfilled exactly. The model developed has been extended to include 
correlations in the quasiparticle pair entering 2q®phonon configuration, i.e. to include two- 
phonon configurations. Finally, the method of elimination of the spurious states within the 
self-consistent QTBA has been considered. 

This work was supported by the DFG (Germany) under grant No. 436 RUS 113/806/0-1. 
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